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Abstract. A new calculus based on fractal subsets of the real line is formulated. 
In this calculus, an integral of order a, < a < 1, called _F Q -integral, is denned, 
which is suitable to integrate functions with fractal support F of dimension a. 
Further, a derivative of order a, < a < 1, called ^"-derivative, is defined, which 
enables us to differentiate functions, like the Cantor staircase, "changing" only on 
a fractal set. The F a -derivative is local unlike the classical fractional derivative. 
The .F^-calculus retains much of the simplicity of ordinary calculus. Several 
results including analogues of fundamental theorems of calculus are proved. 

The integral staircase function, which is a generalisation of the functions like 
the Cantor staircase function, plays a key role in this formulation. Further, it 
gives rise to a new definition of dimension, the 7-dimension. 

F a -differential equations are equations involving ^"-derivatives. They can 
be used to model sublinear dynamical systems and fractal time processes, since 
sublinear behaviours are associated with staircase-like functions which occur 
naturally as their solutions. As examples, we discuss a fractal-time diffusion 
equation, and one dimensional motion of a particle undergoing friction in a fractal 
medium. 



1. Introduction 

It is now a well established fact that fractals can model many structures found in 
nature The geometry of fractals is also a well explored subject ^El ED El IHj- 

Fractals are often too irregular to have any smooth differentiable structure defined 
on them, and render the methods and techniques of ordinary calculus powerless or 
inapplicable. For example the derivative of a Lebesgue-Cantor staircase function is 
zero almost everywhere and therefore this function is not a solution of an ordinary 
differential equation. Consequently, ordinary calculus does not equip us to handle 
problems such as fractal time random walks, anomalous diffusion, dynamics on 
fractals, fields of fractally distributed sources etc., by setting up and solving ordinary 
differential equations. 

During recent times, a few approaches have been developed to deal with various 
aspects of the problems mentioned above. 

Several authors have recognized the need to use fractional derivatives and 
integrals to explore the characteristic features of fractal walks, anomalous diffusion, 
transport, etc. by setting up fractional kinetic equations, master equations and so 
on Ul 00 EH El ID] El Q3I ■ Fractional derivatives are nonlocal operators and often are 
suitable for modelling processes with memory but not always suitable to handle the 
local scaling behaviour e. g. the behaviour of fractal functions. In |14U15llTr)lll7j this 
problem was circumvented by renormalising fractional derivatives and constructing 
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local fractional operators. This was further pursued in ^JE]- A particular success 
of this approach was the demonstration of the striking fact that fractal and multifractal 
functions can be differentiated upto an order (fractional) determined by the Holder 
exponent of the function (or dimension of its graph). In particular, Weierstrass' 
nowhere differentiable function was shown |141 1TB] to be differentiable upto order 
(1 — 7), if (1 + 7) is the box dimension of its graph. 

Another remarkable development is analysis on fractals. Many important ideas 
and applications are developed in the realm of analysis on fractals. This approach has 
been extensively used for the treatment of diffusion, heat conduction, waves, etc. on 
fractals — see j2U| El E21 EH an d several references therein. 

There is a further beautiful developement using a measure-theoretical 
approach |24l I25| . It consists of defining derivative as the inverse of the integral 
with respect to a measure and defining other operators using the derivative. This 
avoids the dependence on the structure of the underlying fractal. 

While all these themes have increased our understanding and brought out 
many beautiful connections, a direct and simple approach involving fractional order 
operators on fractal sets is only moderately explored. Even though measure theoretical 
approach is elegant, Riemann integration like procedures have their own place. They 
are more transparent, constructive, and advantageous from algorithmic point of view. 
It indeed seems possible to develope such an appropriate calculus, tuned to these 
requirements. In the present paper, the first of a series devoted to these ideas, 
we undertake a systematic developement of calculus on fractal subsets of real line, 
involving integrals and derivatives of appropriate orders a £ (0, 1]. A brief glance at 
the table in | Appendix C| would reveal that much of the simplicity and intuitive appeal 
of ordinary calculus can be retained. 

Differential equations of the form D a f(x) = g(x), where D a is a local differential 
operator of order a, < a < 1 were considered in ^B]. It was argued that such 
equations can have meaningful solutions if the dimension of the support of g is a. 
Then the fractional integral of the characteristic function of a Cantor set should be 
the corresponding Cantor staircase function, while the local fractional derivative of the 
latter should be the former. This was done using a Riemann integral like prescription. 
In particular a fractional diffusion equation was shown to have subdiffusive solutions. 
This enabled one to identify a new exact solution of Chapman-Kolmogorov equation. 

This paper formulates the ideas of integral and derivative of order a, < a < 1 
based on a (fractal) set F C R, indicated in ^J- We call them F"-integral and 
^"-derivative respectively. 

The organisation of the paper is as follows: We begin in section [3 by defining a 
mass function 7 a (F, a, 6) and integral staircase function Sp of an order a for a set F. 
The mass function j a (F, a, b) gives us the content of a set F in the interval [a, b) C R. 
Its definition is based on Riemann-like sums. The construction can be compared to 
the definition of Hausdorff measure, except that the covers are more restrictive: they 
are in the form of finite subdivisions of [a, b}. Though 7° is not a measure due to this 
simplification, it turns out to be proportional to Hausdorff measure for compact sets. 
The integral staircase function S F (x) — "f a (F, a, x), obtained from the mass function 
by fixing a, is a generalization of the well known functions such as the Lebesgue- 
Cantor staircase (or the Devil's staircase) functions. The definitions of F a -integral 
and F"-derivative use the quantity {Sp(y) — Sp(xj) in place of the length (y — x) of 
the interval [x, y]. In this respect, the definition of F Q -integral is similar to that of 
Riemann-Stieltjes integral |2fil I27j . 
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In section [21 it is shown that the mass function leads to a new definition of 
dimension called 7-dimension, which is finer than the box dimension, though not 
as fine as the Hausdorff dimension. In later sections it is seen that the ^"-integral 
or the ^"-derivative give meaningful results when the 7-dimension of the underlying 
fractal F is the same as a. 

In section ["£] it is shown that the Hausdorff measure and the mass function agree 
for compact sets upto a proportionality constant. Using this property, the staircase 
function is calculated for the middle | Cantor set. 

Several sets can give rise to the same staircase function. A representative set from 
such an equivalence class of sets, with nice properties, needs to be chosen for defining 
the F Q -derivative and proving the analogues of fundamental theorems. Section [5] 
assures the existence and uniqueness of such a set, called an a-perfect set, associated 
with a staircase function. 

We develope the rest of the theory in a way analogous to the standard 
calculus |261 1271 . In section we introduce notations for limit and continuity 
using the topology of F with the metric inherited from R. This is done in order to 
distinguish them from limit and continuity on R. 

The ordinary integral of functions with fractal support F C R is zero or undefined 
depending on the definition of integral (Lebesgue or Riemann) and the nature of the 
support. The ^"-integral defined in section [7J suits the needs of integration of such 
functions. It is further shown that the F Q -integral of the characteristic function of F 
is the staircase function associated with F as indicated in """""] • 

Functions representing intermittent phenomena or fractal time evolution typically 
"change" only on a fractal. The Cantor staircase function is an example. The F a - 
derivative defined in section |""| is best suited to quantify the "rate of change" of such 
functions. This derivative is local unlike fractional derivatives |29l 1301 l3*Tl I32j . It is 
also not any kind of average derivative as in |331 1341 1351 1361 13*7] . It is more like the 
first order derivative in ordinary calculus, which makes its dynamical interpretation 
possible. Further in the same section, it is shown that the F"-derivative of the staircase 
function of F is the characteristic function of F. Analogues of Rolle's theorem, the 
law of the mean and Leibniz rule are discussed. In section [3] we prove the analogues 
of fundamental theorems of calculus. The formula for F"-integration by parts follows 
thereby. 

The definitions of F Q -integral and the F"-dcrivative reduce to those of usual 
Riemann integral and first order derivative respectively when a = 1 and F = R. 

In section IT""* we discuss examples including subdiffusion and motion of a particle 
undergoing friction in a fractal medium, and demonstrate the use of ^"-differential 
equations as their models. 

As an example, the F"-integral of f(x) — xxc( x ) f° r the middle | Cantor set C 
is calculated in Appendix A| Repeated F"-derivatives and ^-integrals are discussed 
in |Appcndix B where we also calculate ^"-derivatives and F"-integrals of powers 
(Sp(x)) n . A few analogies between classical calculus and F Q -calculus are tabulated 
in |Appcndix C| 

We begin by defining the integral staircase function. 

2. The mass function and the integral staircase 

Let F be a subset of the real line. In most of the interesting cases discussed below, F 
would be a fractal. In this section we formulate the notion of the content or a-mass 
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of F in an interval [a, b], i. e. mass of F PI [a, b], of order a, < a < 1. 
In all the following discussion, < a < 1 unless stated otherwise. 

Definition 1 The flag function 0(.F, /) for a set F and a closed interval I is given by 

6 f I / 1 if ^n/^0 

1 otherwise 

Definition 2 j4 subdivision f^ ,6] (or just P) of the interval [a,b], a < b, is a finite 
set of points {a = xo, x\, ■ ■ ■ , x n — b}, Xi < X4+1. Any interval of the form [xi, Xi+i] 
is called a component interval or just a component of the subdivision P. If Q is any 
subdivision of [a, b] and P C Q, then we say that Q is a refinement of P. If a = b, 
then the set {a} is the only subdivision of [a, b] . 

Definition 3 For a set F and a subdivision P[ a ,b] , a < b, 

° a l F > p ] = E u XT) ^ [^^m])- (2) 

f — ' 1 (a + 1 



i=0 



If a = b, we define o~ a [F,P] to be zero. 

We note that the sum in © contains a contribution from a component interval if and 
only if that component contains at least one point of F. Further, cr a [F, P] > for any 
set F and subdivision P of [a, b]. 

We remark that this definition, and in particular the factor 1/T(a + 1) and the 
use of finite subdivisions, has been motivated by local fractional calculus |lf)U17| . 

Now we introduce the coarse-grained mass: 

Definition 4 Given 5 > and a < b, the coarse-grained mass "f$(F, a, b) of F D [a, b] 

is given by 

1 f{F,a,b)= inf ■ a a [F,P] (3) 

{"[a,6] : l"l<0} 

where 

\P\ = max (x i+1 - Xi) (4) 

0<i<n-l 

for a subdivision P , and the infimum in is taken over all subdivisions P of [a, b] 
satisfying \P\ < S. 

Eventually, a limit of 7^(^,0,6) as 5 — > is taken in definition [S] below. But 
before that we examine some important properties of 7^(^,0,6). 

Let a < b and Si < 62- Then 7^ (F, a, b) is the infimum of cr a [F, P] over a smaller 
class of subdivisions than ^(F, a, b). Thus: 

Lemma 5 Let a < b and Si < S2 ■ Then 7^ (F, a,b) > 7!? (F, a,b). 

The following lemma shows that 7^ (F, a, b) is a monotonic increasing function of 
b and a monotonic decreasing function of a. (Throughout the paper we distinguish 
between monotonic and strictly monotonic.) 

Lemma 6 Let 5 > and a < b < c. Then, ^f{F,a,b) < ^f{F,a,c) and 
Tf{F,b,c)<>tf(F,a,c). 
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Proof: Let e > 0. Then according to the definition 0] of r )f{F 1 a, c), there exists a 
subdivision P\ a , c ] — { x o — 

= c} such that \P\ < S and 

a a [F 1 P]< 1 £(F,a,c)+e. 

Let Q[ a ,b] = {x e P : x < b} U {b} i. e. Q[ a ,b] = {Vo, yi, ■ ■ ■ , y m } where y t = x, if x, < b 
and y m = b. 

It follows that \Q[ a ,b] \ < \P[a,c]\ < <5 and 0(F, [y m -i,y m ]) < 9{F, [x m -i,x m ]) since 
[ym-i,ym] C [x m -i,x m ]. Therefore, 

o~ a [F, Q [aM ] < a a [F, P M ] < ^{F, a, c) + e. 
But "fg(F, a, 6) < o- a [F, Q] and e is arbitrary, hence 

-ti(F,a,b) < 7 f(F,a,c). 
which completes the proof of the first part. The second part follows in a similar way. 

Theorem 7 jf(F,a,b) is continuous in b and a. 

Proof: We prove the continuity of jg(F, a, b) in b (with S, a and a fixed). Continuity 
in a can be proved in a similar manner. 
Given e > 0, let 

A' = (eT(a + 1))= and A = min(A', S). 

For ei > 0, there exists a subdivision P, such that |P| < 6 and 

a a [F,P] <^{F,a,b) + ex. 

Now Q = P U {6 + A} is a subdivision of [a, 6 + A]. Therefore, 

7 5 Q (F,a,& + A) <a a [F,Q] 

^a"[F,P]+g(F,[b,b + A]) A " 

1 (a + 1) 

< a Q [P,P] +e 
<7, Q (P,a,6)+ ei + e . 

As ei is arbitrary, we get "fg(F, a, b + A) < j£(F,a,b) + e. As jg(F, a, b) is a 
nondecreasing function of b, 

tf(F,a,b + t)<-y£(F,a,b)+e 

for < t < A. 

Summarizing, given e > 0, there exists a A > such that 

c - 6 < A => j£(F, a, c) - 7 , Q (P, a, b) < e 

which implies that J^{F, a, b) is continuous in b from right. The continuity from left 
follows on the replacement of b by b — A and of b + A by b in the above proof. • 

As mentioned earlier, the mass function is the limit of the coarse-grained mass as 

Definition 8 The mass function 7 a (P, a,b) is given by 
j a (F,a,b) = ^m^(F ) a,b). 
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We note that since lg(F, a, b) increases as S decreases, 7 Q (F, a, b) always exists and is 
a non-negative number, which may possibly be +00. 

Another simple observation is that if F n [a, b] = 0, then 7 ^(F, a, b) — for any 
5 > 0, and consequently 7 Q (F, a, 6) = 0. This result can be extended so that it also 
applies to an open interval (a, 6): 

Lemma 9 If F n (a, b) = 0, then a, 6) = 0. 

Proof: If F n [a, 6] = 0, then the result is obvious. If not, then F n [a, 6] contains one 
or both of a and b. In that case, given e > 0, we can choose a subdivision P of [a, 6] 
such that 



(#1 - xq) and (x n - x n -\) < 



er(a+l) 
2 



where {xq, . . . , x n } are points of P. Then, o~ a [F, P] < e since [xi, x„_i] (~1 P = 0. But 
as e is arbitrary, it follows that 7"(P, a, 6) = for any <5 > 0, so that 7 Q (P, a, 6) = 0. 



A property, desired of a mass function, is additivity. The following theorem 
asserts this. 

Theorem 10 Let a < b < c and 7 a (P, a, c) < 00. Then 

7 "(P, a, c) - 7 Q (P, a, b) + 7 Q (P, 6, c). (5) 

Proof: Given <5 > 0, let Pi be any subdivision of [a, b] and P2 be any subdivision 
of [b,c] such that |Pi| < 5 and | P-2 1 < 5. Then, Pi U P2 is a subdivision of [a, c], 
I Pi UP 2 | < <5, and 

ct" [P, Pi U P 2 ] = ct" [P, Pi] + ct q [P, P 2 ] . 

Taking infimum over all subdivisions Pi and P 2 such that |Pi| < S and |P 2 | <S, and 
noting that not all the subdivisions of [a, c] can be written in the form Pi U P 2 , where 
Pi is a subdivision of [a, b] and P 2 is that of [b, c] , we get 

ls{F,a,c)< inf <7 a [P,PiUP 2 ] 
|Pi|<«,|P 2 |<« 

= 7 f(F,a,b)+ 7 f(F,b,c). (6) 

Now for every subdivision P[ a , c ], < <^ we can construct a subdivision 
P' = P U {b}. Obviously \P'\ < 5, and P' = P 1 U P 2 where Pi is a subdivision 
of [a, 6] and P 2 is a subdivision of \b, c\. 

Let P = {x ,x u ...,x n }. lib G P, then P = P' and cr"[P,P] = ct q [P,P']. 
Otherwise, let [xfc,Xfc+i] be the interval which contains &. Thus, 

a«[p,pu {&}] - <7«[p,p] = e(F, M) ^~ x + k) ^ +e(F, [MwD^Y^f 

-0(P,[x fc ,x fe+ i])- (Xfc+1 " Xfc)a 



r(a + l) 
Hence, 

^[F.PU^l-^lF.P]^^. 
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This implies that 

a a [F,P} + >a a [F,PU{b}} 
T(a + 1) 

= a a [F,Px]+u a [F,P 2 ] 
>tf(F,a,b)+>tf(F,b ) c) 
for all P. Thus if we take infimum over all subdivisions P such that \P\ < <5, we get 

Is (F, a, c) + r( — t > 7? OF, a, 6) + 7? 6, c) (7) 
1 (a + 1) 

From 10 and J7J) and taking the limit as 5 — ► 0, we get the result. • 

Since each term in is nonnegative for a < b < c, an immediate consequence is 

Corollary 11 j a (F,a,b) is increasing in b and decreasing in a. 

The next theorem states that j a (F,a,x) takes all values in the range 
(0, 7 a (F,a,6)) for x £ (a, b). 

Theorem 12 Let a < b and let j a (F,a,b) / be finite. Let y be such that 
< y < 7 a (-F, a, b). Then there exists c, a < c < b, such that "f a (F, a, c) = y. 

Proof: Let z = 7 Q (F, a, 6) — y. 

Given a (5 > 0, consider the set of all points x of [a, b] such that 7 J*(-F, x, 6) < z. 
This set is an interval of the form [sg, b] for some sg, a < s$ < b, because r )f{F 1 x, b) is 
continuous (theorem[7|) and decreasing in x fcorollarv lllfl . Since j"(F,x,b) increases 
as 6 decreases (lemma |5J|, ss increases as 5 decreases. 

Similarly the set of all points x of [a, b] such that jf(F, a,x) < y is an interval of 
the form [a, t^], a < t$ < b, and tg decreases as S decreases. 

Let x £ (a, b). Then by theorem 1 101 

7 a (F,a 1 b)= 1 a (F,a,x)+ 1 a (F,x,b)>^(F,a,x)+ 1 nF,x,b). (8) 

As y, z < 7 Q (i 7 ', a, b), there exists a Sq > such that 5 < Sq implies that 
jf(F, a, b) > y, z. In the rest of this proof, we only consider S < Sq without mentioning. 

Since jg(F, a,b) > y and j"(F, a, u) is continuous and increasing in u, there exists 
an x £ (a, b) such that "fg(F, a, x) = y. This implies that x £ [a, tg]. Further, from 
iJSJ, it follows that 

z = 7 Q (F, a, 6) - y = j a (F, a, b) - ^(F, a, x) > 7 f(F, x, b) 

implying that x also belongs to [sg, b]. This can happen only when sg <tg. 
Thus for each 5 there exists an interval [sg , tg] such that 

x £ [sg, tg] 7™ (F, x,b) < z and jg(F, a, x) < y. 

Let s = snp 0<s< g sg and let t = info<5<5 tg. Now sa increases and tg decreases as 6 
goes to zero, but as sg < tg for any 5. Thus s < t and 

[m] = n is*,**]- 

0<<5<<5 

Consequently x £ [s,t] implies jg(F,x, b) < z and J?(F, a, x) < y for any i5. Hence 

x £ [s, t] j a (F, x,b) <z and j a {F, a,x) <y (9) 
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But as 7 a (i 7 ', a, x) + 7 a (-F, x, b) = j a (F, a,b) — y + z, the inequalities in @ must be 
equalities. Thus for a given y, < y < "f a {F, a, 6), there exists a set [s, t] C [a, b] such 
that a; G [s, t] l a {F, a,x) = y which completes the proof. • 



Corollary 13 If j a (F, a,b) is finite, ^"(F, a, x) is continuous for x G (a, 6). 

This can be proved using the monotonicity of 7 Q (-F, a, b) in a and b. 

Remark: The implication of this result is that no single point has a nonzero mass, or 

in other words, the mass function is atomless. 

The scaling and translation properties of the mass function are similar to those 
of Hausdorff measure: 

Theorem 14 For FcR and A G R, let F + X denote the set 

F + \ = {x + \:xeF} 
and let XF denote the set 

XF = {Ax : x G F}. 

Then, 

(i) Translation: 

7 Q (F + X,a + X,b + A) = j a (F,a,b) 

(ii) Scaling (X > 0): 

j a (XF,Xa,Xb) = X a j a (F,a,b) 

Remark: If the set F is self-similar so that XqF n [Aoo, Aofe] = F n [Xqo, Xob] for a 
particular Ao, then the scaling property can be rewritten as 

7 Q ( J F 1 ,A a,A 6) = A^7 Q (^a^)- 
An example is the middle | Cantor set C, with a = 0, b = 1 and Ao = ^r- 

Now we introduce one of the central notions of this paper, viz. the integral 
staircase function for a set F of the order a. This function, which is a generalization 
of functions like the Lebesgue-Cantor staircase function, describes how the mass of 
F n [a, b] increases as b increases. 

Definition 15 Let ao be an arbitrary but fixed real number. The integral staircase 
function Sp(x) of order a for a set F is given by 

, ^ \ l a (F,ao,x) if x > ao 

S»(x) = < „ , ' . ~ (10) 

I —7 (F,x,ao) otherwise. 

The number ao can be chosen according to convenience. A few properties of Sp(x) 
which are restatements of the corresponding properties of the mass function j a (F, a, b) 
are as follows. 

Theorem 16 Let F be a subset o/R 7 and let < a < 1. If "f a (F, a,b) is finite, then 
for all x, y G (a, 6) such that x < y, the following statements hold: 

(i) Sp{x) is increasing in x. 

(ii) If F (~l (x,y) — 0, then Sp is a constant in [x,y]. 
(tit) Sp-(y)-Sp-(x)=i a (F,x,y). 

(iv) Sp is continuous on (a, 6). 

As an example, we calculate and show the graph of Sq for the middle ^ Cantor 
set C in the section 0] after discussing some results required to calculate it. 
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3. The 7-dimension 

We now consider the sets F for which the mass function j a (F,a,b) gives the most 
useful information. Due to the similarity of the definitions of the mass function and 
the Hausdorff outer measure 0] |5J E] , one might expect that the mass function 
can be used to define a fractal dimension. It is indeed the case. If < a < f3 < 1, 

AF,P]<\Pf- a ^[F,P]^±^ 

so that 

^ { F,a,b)<S^(F,a,b)^±^ 

Thus in the limit as d — * 0, we get 

7 /3 (F, a, b) = provided j a (F, a, b) < oo and a < (3. 

It follows that j a (F, a, b) is infinite upto certain value of a, say cto, and jumps down 
to zero a > cto (if ao < 1). We call this number the 7-dimension of F. j a ° (F, a,b) 
itself may be zero, nonzero finite, or infinite. To make the notion of dimension precise, 

Definition 17 The 7- dimension of F n [a, b], denoted by dim 7 (i 7 ' n [a, b]), is 
dim 7 (F n [a, b}) = inf{a : j a (F, a, b) = 0} 
= sup{a : ^{F, a, b) = 00} 

Now we compare the 7-dimension with the Hausdorff dimension and the box 
dimension. As the definition of Hausdorff measure involves arbitrary countable covers, 
it is expected that the Hausdorff dimension be finer than 7-dimension. This is shown 
to be the case below: 

Let TL^(E) denote the coarse grained Hausdorff measure of a subset E of R, 
and T-l a (E) denote the Hausdorff measure. Let P be a subdivision with \P\ < S. 
Those components [x^Xi+i] of P for which 9{F, [xi,Xi+\]) is nonzero, form a 5-cover 
of Fn [a,b]. Thus, 

<J a [F,P}>— ^— Hf(Fn[a,b}). 
T{a + 1) 

Since this is true for any P such that \P\ < 5, it follows that 

if{F, a, b) > —J—Hf(F n [a, b]) 
L(a + 1) 

for each S > 0. So taking limit as 8 — > 0, 

l a {F,a,b) > —±— H a (Fn[a,b}). (11) 

which also implies 

dim w (F n [a, b} ) < dim 7 (F n [a, b} ). 

There exist sets for which the two definitions give different results. For example, 
if Q denotes the set of rational numbers, then dim-^(Q n [0, 1]) = 0, while dim 7 (Q n 
[0, 1]) = 1. However, it will be shown in the next section that the two dimensions are 
equal for compact sets. 
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Next we compare the 7-dimension with the box dimension. Let dim 7 (i 7 'n [a, b]) = 
a. Then ^(F, a, 6) diverges for any (3 < a. Thus for any k > 0, there exists <5o > 
such that S < do => Jg(F, a, b) > k. 

Let P be any subdivision such that \P\ < 5, and let N$(F n [a, 6]) be the number 
of nonzero terms in the sum o~ a [F, P] . Then, for arbitrary but fixed k > and S < 6q , 

fc<7f(F,«,b)< fC9 + l) 
where < j3 < a < 1. Thus, 

ln(fc) < lnJV 5 (F n [a, 6]) + /31n(£) - ln(r(/3 + 1)) 

which implies 

-/31n(<5) < In Ns(F n [a, 6]) - In(jfc) - ln(r(/3 + 1)). 
Dividing by — ln(<5) (which is positive for 5 < 1), 

ln(iNfr(F n [a, fe])) ln(fc) - + 1)) 

P " - ln(5) - ln((5) 

Taking limit as S — > and noting that the first term is the definition of the box 
dimension dims(i 7 ' n [a, b]) in the limit and the denominator of the second diverges, 
we get 

n ,. r ,in ln(iV 5 (Fn [a, 6])) 

(3 < dim B Fn o, ft ) = lim v ' v - .. L . J;; . 

5^0 — ln(o) 

This is true for any /3 < a = dim 7 (i 7 ' n [a, 6]), so that 

dim 7 (F n [a, 6]) < dim B (F n [a, 6]). 

As an example in which box dimension and 7-dimension differ, we consider the 
set K — {0, 1, 3, 5, • ■ ■}• It is known 0] that dims K — 0.5, and dim^ K = 0. Further, 
as K is compact, dim 7 K = dim-^ K according to the corollary ^j] to be proved in the 
next section. 

Thus the 7-dimension is finer than the box dimension, but not than Hausdorff 
dimension. Specifically, the 7-dimension is unaffected by clusters of points unlike the 
box dimension, as the above example of the set K demonstrates. On the other hand, 
it is sensitive to countable but dense sets such as rationals. 

The modified box dimension 4 involves countable subsets of the set under 
consideration, to circumvent the above mentioned problems in box dimension. 
Therefore the former is finer than 7-dimension. But the 7-dimension is simpler to 
calculate. For similar reasons, the 7-dimension is simpler than packing dimension 4 . 

Summarizing, if E C R, then 

dims E > dim 7 E > dim^/s E > dim-^ E. 



4. The mass function and the Hausdorff measure 



For compact sets, any open covers, and in particular countable open covers, can be 
replaced by finite subcovers. So the mass function and Hausdorff measure arc expected 
to be proportional. Theorem 1181 below shows that this is indeed the case. In what 
follows, Hg(F) denotes the coarse grained Hausdorff measure _3] of a set F C R, and 
TL a (F) denotes the Hausdorff measure of F, both of order a. 
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Theorem 18 For a compact set F C R, 

7tt(f - n '^ r( a + i) r(fn|a ' tl) - 

Proof: For 5 > 0, let {Ai, i = 1, 2, . . .} be any countable cover of F n [a, b] such that 
diamAi < I for all i. The sets Ai need not be open or closed. Then 

Hs /2 (F n [a, 6]) < ^(diamAA Q . 

i 

Consider closed intervals Bi = [ui,Vi] where ut — inf Ai and Vi = sup Ai. Then 
A.; C Si and diami?i = diamAi. Thus {Bi} forms a cover of F R [a, 6] and 

^(diam^)" = ^(diamA^ > H? /a (F H [a, 6]). 

z i 

Given e e (0, (5/2)"), let {CJ^ be the open intervals 
Ci = 

The class {Ci}^ thus forms an open cover of F D [a, b] and 
diamCi = diamAi + ° < ^ 

so that 

^(diamar = E (diamA, + (^) ' 

A simple consequence of Jensen's inequality |SH1 > which for the case of two variables 
assures that (si + S2Y < s\ + s| for s±, s 2 > and < t < 1, is that 

^(diamC,) < E( diamA *) Q + e X^ ^ E( diamA *) Q + e ' ( 12 ) 

We now show that a finite cover consisting of closed intervals can be constructed. 
As F is compact, so is F fl [a, b}. Thus a finite subset of {d} covers F fl [a, &]. We 
denote this finite subcover by {Z?i,i = 1, . . . , n}. The Di are open intervals of the 
form (cij, 6j). Without loss of generality we can choose this finite subcover {Di} such 
that Di <f_ Dj whenever i ^ j. Further, the sets are labeled such that <ij < a i+1 . But 
as Di <f_ -Dj+i and D i+ i <f_ Di, it implies that ctj < a J+ i and bi < bi + \. 

Now we consider the closures Di of Di. As {-Di} is a finite subcover out of {Ci} 
and have the same diameters as Di, it follows from (|12|) that 

^(dianxDA* < ^(diamA,)" + e. 

Let 7i = Z?i and ii = £)j\Z)j_i for 2 < i < n. The collection {ii} forms a finite cover 
of F fl [a, 6] by closed intervals, and 

^(diam/i)" < ^(diamA,)" + e. 

The closed intervals ij share at the most endpoints. The set of all the endpoints of Ii, 
1 < i < n forms a subdivision P of [a, b] which can be refined to a subdivision Q such 
that \Q\ < S and 

r(a + 1)ct q [F, Q] = ^(diam/i) Q < ^(diamAi)" + e. 
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Therefore, 

T(a + 1)7?(F, a, b) < ^(diam^)" + e. 

Since this relation holds for any countable cover {Ai} of FU [a, b] such that diamAi < 
8/2 and for arbitrary e > 0, it follows that 

T(a + lh^F, a, b) < H% /2 (F n [a, b}) < H a (F n [a, b}). 
Consequently in the limit as 5 — > 0, 

r(a + l)7 a (F, a, 6) < n [a, 6]). (13) 

Equations l|ll|) and l|13[l together imply the required equality. • 

Corollary 19 If F C R is compact, then dim 7 i 7, = dim^ i<\ 

Example: We now discuss an important prototype example of S^x). Consider 
the middle 1/3 Cantor set C (hereafter referred to as the Cantor set). This set is 
compact and has a Hausdorff dimension a — log(2)/ log(3). Thus by theorem 1181 
U a {C n [a, b]) = T{a + l)"f a {C,a,b) and dim 7 C = dim H C = a = log(2)/log(3). 
Using the self-similarity of C and the monotonicity as well as scaling and translation 
properties of the mass function (theorem [2J|, we can calculate Sp at each point. A 
graph of T(a + 1)Sq{x) i. e. T(a + 1)7" (C, 0, x) against x is shown in figure Q This 
is the Lebesgue-Cantor Staircase function. 



- 

0.9 - I - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 1. T(a + 1)Sq{x): The integral staircase function for the Cantor set 



5. a-Perfect sets 

The correspondence between sets F and their staircase functions Sp is many to one. 
For example, Sq, = Sq where C is the middle | Cantor set and C = C\{|}. 
Intuitively, it can be said that as the mass function is atomless, removing a single 
point from C does not change its value. 

Another example is a set D — C U E where E C (|, | ) satisfying dim 7 E < a, 
a = ln(2)/ln(3). Then it can be seen that Sjj(x) — Sq(x) for all x. Thus adding a 
lower dimensional set need not change the value of the staircase function either. We 
call the sets giving rise to the same staircase function as staircasewise congruent: 
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Definition 20 Let FcR and G C R be such that dim 7 F — dim 7 G — a. a e (0, 1]. 

T/ien _F and G are said to be staircasewise congruent if S F (x) and Sq(x) are finite 
and equal for all x£R. 

This congruence being an equivalence relation, we denote the equivalence class of 
sets containing F by Sp. Thus if G is in S F , then Sq = Sp and Sq = £p. 

The above examples intuitively suggest that not all points or subsets contribute 
to the staircase function. Now we proceed to select a representative set out of the 
equivalence class which, intuitively speaking, has exactly those points at which Sp 
"changes". To choose only the points where a function "changes", we need the 
following definition. 

Definition 21 We say that a point x is a point of change of a function f , if f is not 
constant over any open interval (c,d) containing x. The set of all points of change of 
f is called the set of change of f and is denoted by Sch/. 

Thus Sch/i = if /i = constant, while Sch/2 = R if fz(x) — x. More importantly 
we note that if G G then S*g(x) = Sf{x) and therefore Sch(Sg) = Sch(S£). In 
other words, Sch(S^) is determined by the equivalence class Sp. 

The following theorem states that Sch(S^) itself belongs to Sp. 

Theorem 22 Let F C R be such that Sp(x) is finite for all x G R for a — dim 7 F 
and H = Sch(S^). Then H belongs to £f i. e. S% = S F . 

Proof: It will be sufficient, in view of definition 1151 to prove that j a (H,a,b) = 
7 Q (F, a, b) for any a, b G R. 

We begin by noting that for u < v, if F H [u,v] = 0, then -f a (F,u,v) = 0. 
Consequently Sp is constant on (u,v) implying (u,v) PI H = 0. Then for any e' > 
such that u + e' < v — e', 

6{H, [u + e',v-e']) = 0. (14) 

Next, let 5 > 0. Then given e > 0, there is a subdivision P[ a ,b] — {no, Ml-, • ■ ■ 1 Un} 
such that \P\ < S and 

o- a [F,P] < 7 «(F, 0> 6) + 1. (15) 

If 9(F, I) — I for all components L of P, then certainly 

a a [H,P]<a a [F,P]. (16) 
Otherwise, let K be the set of all points of the form 



2n J ' \ 2n J 



(17) 



where c and d are the endpoints of those components I of P such that 0(F, I) = 0. 
Then Quu — PUK is a refined subdivision and \Q\ < \P\. If I = [c, d] is a component 
of P such that 9(F, L) = 0, then it contains three components of Q, viz. [c, c'], [c 1 , d'\ 
and [d' , d) where d and d! are given by (TTJ. The term in cr a [H, Q] corresponding to 
[c',g?'] is zero according to 1)14(1. and the remaining two contribute at the most e/2n 
each. If 9{F , I) ^ 0, then I is also a component of Q and the term corresponding to / 
in a a [H, Q] is either zero or is exactly the same as the corresponding term in a a [F, P] . 
Therefore, 

a a [H,Q] <a a [F,P]+l (18) 



Calculus on fractal subsets of real line - I: formulation 14 

Thus from JTHJl, l(Po|) and (jT%)l . we see that there exists a subdivision Q, such that 
\Q\ < 5 and 

<J a [H,Q] <tf(F,a,b) + e. 
As e is arbitrary, we see that 

7 ?(^«^)<7 5 Q (^«,6) (19) 
Now we wish to rule out the possibility that 

1 ?(H,a,b)< 7 f(F,a,b). (20) 

Suppose that (|2(J[1 is true. Then there exists a subdivision Pi = {xq, . . . , x n } such 
that | Pi | < S and 

a a [H,Pi] < -ff(F,a,b). (21) 

From 10 we know that 

n-l 

7? <$>i(*>i,*i+i) 

i=0 

so that 

n-l 

a Q [ff,Pi] < ^ 74(^,^,^+1 )■ 

i=0 

For this equation to hold, there must be at least one k, < k < n — 1, such that 

fffc+l Zfc) g [a; feja ; fc+1 ]) < 7^(F,x fc ,2: fc+ i) < 7 Q (F,a; fc ,2:fc + i). 
1 (a + 1) 

As every quantity is non- negative, it follows that 7 a (F, x kl Xk+i) > 0, and Sp is not 
constant in [x kl x k +{\. Therefore Sch(5J?) n [xfc,Xfc+i] ^ i. e. H n [xfc, Xfc+i] ^ 0. So 
6>(if, [x fc ,a;fc + i]) = 1 and 

(xk+i - Xk) 



< l&{ F ,%k,Xk+l)- 

on of [x k ,x k+ i] t 
a a [F,Q]<^(F,x k ,x k+1 ) 



r(a + i) 

As Q = {xfc, Xfc + i} is a subdivision of [x k , x k +i\ such that \Q\ < 5, 

{x k+ i - x k ) a 



r(a + i) 

which is a contradiction by the definition of 7"(F, Xfc, Xfc+i) implying that our 

assumption (121 H is wrong. Thus, (|19|l is an equality for any S > 0, and therefore 
7 Q (ff,a,6) = 7 Q (P,a,6). • 



Lemma 23 Let F C R be such that Sp(x) is finite for all x G R /or a = dim 7 F. 
TTien t/ie set if = Sch(S'^) is perfect i. e. H is closed and every point of H is its limit 
point. 

Proof: Let y be a limit point of H . Then any open interval (c, d) containing y contains 
a point z of H = Sch(S^). Therefore Sp is not constant on (c, d). Hence y E H 
implying that F is closed. 

If x £ H is not a limit point of H, then there exists an open interval (c, d) 
containing x but no other point of H so that F n (c, x) — and F l~l (x,d) = 0. 
This implies that Sf(d) - Sf(c) = (Sf(d) - Sf(x)) + (Sf(x) - Sf(c)) = due to 
theorem lltjf ii). Therefore x is not in H = Sch(Sp) which is a contradiction. • 

Now we choose Sch(S'^) as the "canonical" representative of £p\ 
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Definition 24 Let F C R be such that Sp(x) is finite for all x E R for a = 
dim 7 F. Then the set Sch(Sp) is said to be a-perfect, and is said to be the a-perfect 
representative of £p. 

Thus, every £p contains a unique a-perfect set. The next theorem states that it is 
the minimal closed set in £p. 

Theorem 25 An a-perfect set F is the intersection of all the closed sets G in £p. In 
other words, it is the minimal closed set in £p. 

Proof: Let Q be the class of all closed sets in £p. 
As F is perfect, it is closed. Therefore, 

F D p| G. (22) 
Gee 

Let Go E G and x ^ Go- Then there is an open interval (c, d) containing x but 
no point of G , as G is closed. This implies that Sq (c) = Sq (d) by theorem HTfl nl. 
and further S F (c) — S F (d) as G 6 £ F . Since F is a-perfect, we have F tl (c, d) = 
implying that x £ F. Therefore, x £ F => x e Go for all Go £ G, so that 

F c p| G. (23) 
Gee 

The proof is completed in view of (|22|l and 1|23|) . • 

The following lemma is a useful restatement of the definition of an a-perfect set. 



Lemma 26 Let F C R 6e a-perfect and x £ F. If y < x < z, then either 
S F (y) < Sf(x) or S F {x) < Sf(z) (or both). 

Thus for an a-perfect set F, the lemma assures that if x £ F, then the values of S F (y) 
must be different from Sp(x) at all points y on at least one side of x. 

As an example, we now show that the middle | Cantor set G is a-perfect, for 
a = log(2)/ log(3). Let x £ C. As G is closed, there is at least one open interval (c, d) 
containing x, such that G fl (c, d) = 0. Therefore, 5g is constant on (c, d) implying 
that x i Sch(S'g). Hence, 

Sch(Sg) c G. (24) 

Let a; be a point of C . Then x can be represented by Q.X1X2X3 . . . where Xi 
is the ith digit in the ternary representation of x. As x E C, Xi = or 2. Let 
(c, d) be any open interval containing x. Then there is an integer n > such that 
(x — 3~ 7 \x + 3~") C (c, 0?). Let D be the set of numbers y = 0. 2/12/22/3 • • ■ satisfying 

{Xi l<i<n + l 

or 2 i > n + 1 

where yi is the ith digit in the ternary representation of y. Then D C G is a scaled 
down copy of G by a factor 3~" _1 and D C (x — 3~", x + 3~") C (c, d). Thus using 
the scaling property of j a f theorem HI) !. 

5*g(d) - S*g(c) > — -(3-"- 1 ) = — -2-"- 1 > 

CK ' CK ! ~ T{a + \y ' T(a + 1) 

where a = ln(2)/ln(3). This implies that Sq is not constant on (c, d). Thus, 

x E C => x e Sch(S'g) i. e. 

G c Sch(S'g). (25) 

From (|24|l and (|25|l we see that G = Sch(5g) implying that C is a-perfect for 

a = log(2)/log(3). 
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6. F- continuity 

In this section we introduce the notation for limit and continuity using topology of 
F C R with the metric inherited from R. Our purpose in doing so is to distinguish 
between these notions and ones on R when they both appear together. 

Definition 27 Let FcR, / : R — > R and x G F. A number I is said to be the limit 
of f through the points of F, or simply F -limit, as y — > x, if given any e > 0, there 
exists 5 > such that 

y G F and\y - x\ < S =^ \ f(y) - l\ < e. 

If such a number exists, then it is denoted by 

e = F-limf(y) 

y^x 

This definition does not involve values of the function at y if y £ F. Also, F-limit is 
not defined at points x ^ F. 

We now introduce the notion of F-continuity which is continuity as far as the 
values of the function only on the set F are concerned. 

Definition 28 A function f : R — > R is said to be F -continuous at x G F if 
f(x)=F-\imf(y) 

y^x 

We note that the notion of .F-continuity is not defined at x ^ F. 

It is clear that continuity of / : R — > R at a; e F implies -F-continuity at x. But 
the converse is not true. We consider a few examples: Let C be the middle | Cantor 
set. Then the functions fi(x) — 1 and /^(x) = x are continuous on [0, 1]; they are also 
C-continuous on C H [0, 1]. In contrast, consider f^{x) — xci x ) an d fA{%) — % • Xc( x ) 
where xc{x) is the characteristic function of C. These functions are C-continuous, 
but not continuous. 

Now we define an analogue of uniform continuity: 

Definition 29 A function f : R — > R is said to be uniformly F '-continuous on E C F 
if for any e > there exists 5 > such that 

x £ F,y e E and \y — x\ < 5 \ f(y) - f(x)\ < e. 

It is clear that uniform ^-continuity on E implies _F-continuity on E. The converse 
is true only in certain cases: 

Theorem 30 // a function f : R — > R is F -continuous on a compact set E C F , 
then it is uniformly F -continuous on E . 

7. ^-Integration 

In the definition of i^ a -integral below, values of the function only at the points of 
F are considered. Further, instead of the lengths of subintervals, we consider the 
difference between the values of the staircase function Sp at the endpoints. In this 
respect, F Q -integral is similar to Riemann-Stiltjes integral |261 127j . 

Definition 31 The class of functions f : R > R which are bounded on F is denoted 
by B(F). In other words, 

f G B(F) -co < inf f(x) < sup f(x) < +oo 

X ^ F xGF 
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As the first step, we now define upper and lower sums which approximate the 
value of the F Q -integral. 

Definition 32 Let f e B(F). Let L be a closed interval. Then, 

M[f,F,I}= sup /(*) ifFDl^Q) 
xeFni 

= otherwise 



and similarly 



m[f,F,I}= Jnf^/Oc) 

= otherwise 



Definition 33 Let Sp(x) be finite for x € [a, b]. Let P be a subdivision of [a, b] with 
points Xo, ■ ■ ■ ,x n . The upper F a -sum and the lower F a -sum for the function f over 
the subdivision P are given respectively by 

n-l 

U a [f,F,P] =J2 M [f^ F Axi,x i+1 ]}(S^(x i+1 ) - S^( Xi )) (26) 

and 

n-l 

L a [f,F,P] =Y,m[f,F,[x i ,x i+1 ]](S^{x i+1 )-S^(x i )) (27) 

i=0 

We emphasize the appearance of intersection F n / in the definition of M and m, and 
also the use of (Sp(x i+ i) — S F (xi)) as in a Riemann-Stieltjes sum instead of (x i+ i —Xi). 
From the definition it is clear that 

U a [f,F,P] >L a [f,F,P}. (28) 

The following lemma asserts that with refinements, the upper F a -sum decreases 
and the lower F"-sum increases, both monotonically. 

Lemma 34 Let FcR and f G B(F). If Q is a refinement of a subdivision P, then 
U a [/, F, Q] < U a [/, F, P] and L a [/, F,Q] > L a [/, F, P] . 

Proof: To start with, let P = {x , xi, . . . ,x n } and Q = PL){x'} where x' € (xi, Xi+i). 
Let / = [xi,Xi + i], I' = [xi,x'], and I" = [x',x i+ i]. If there are no points of F in /, 
then M[f, F, I] = M[f, F, I'] = M[f, F, I"] = 0. Otherwise there are two possibilities; 
cither both I' and I" have points of F, or only one of them, say I' without loss of 
generality, has points of F. 
In the first case, 

M[f, F, I'} < M[f, F, I] and M[f, F, I"} < M[f, F, I}. 

Thus we have, 

M[f, F, I](Sp(x i+1 ) - Sp(xi)) 

= M[f,F,I]{{S F {x') - S a F { Xl )) + (S$(x i+1 ) - S^(x'))} 

> M[f,F,I'](S%(x') - S^( Xi ))+M[f,F,I"](S^(x i+1 ) - Sftf)). (29) 

In the second case, only /' has the points of F. Consequently Sp(x i+ i) = Sp(x') 
and M[f, F, I'] = M[f, F, I}. Thus 

M[f,F,I]{S%{x i+1 )-S%{x i )) = M[f,F,I>]{S%{x')-S%{x i )) (30) 
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since M[f,F,I"] = 0. Combining {22J and QQJ, we have 

> M[/,F, /'^(a/) - + M[/,F,/"](^(^+i) - 

Thus 

U a [f,F,Q]<U a [f,F,P}. 

This conclusion can easily be extended for any refinement of P. 
By a similar argument, we can prove that 

L a [f,F,Q] >L a [f, F,P] 

which completes the proof. • 

Lemma 35 If P and Q are any two subdivisions of [a,b], then 
U a [f,FP]>L a [f,F,Q] 

Proof: As P U Q is a refinement of both P and Q, it follows from the above lemma 
and (H that 

U a [f, F P] > U a {f, F,PUQ]> L a [f, F,PUQ}> L a [f, F, Q]. 



We are now ready to define the F Q -integral. 

Definition 36 Let F be such that Sp is finite on [a,b]. For f £ B(F), the lower 
F a -integral is given by 

[ f(x)d F x= sup L a [f,F,P] (31) 

and the upper F a -integral is given by 

[ b f(x)d F x=MU a lf,F,P] (32) 

Both the supremum and infimum are taken over all the subdivisions P of [a,b]. 

The d F x appearing in (j31(l and Ij32(l has no separate meaning; it is just the notation. 
It is obvious that 

rb ~ Jb 

f(x)d a F x< / f(x)d a F x (33) 



Definition 37 If f G B(F), we say that f is F a -integrable on [a,b] if 
f(x)d F x= / f(x)d F x 



In that case the F a -integral of f on [a,b], denoted by J^f(x)d F x is given by the 
common value. 

For future use we note the following obvious and useful criterion for proving 
F Q -integrability: 
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Lemma 38 Let f G B{F). Then f is F a -integrable on [a,b] if and only if for any 
e > 0, there exists a subdivision P of [a, b] such that 

U a [f,F,P]<L a [f,F,P]+e. 

Now we state a sufficient condition for i^-integrability. The sufficient and 
necessary conditions will be discussed in a companion paper. 

Theorem 39 Let F be such that F n [a, b] is compact and Sp is finite on [a, b]. Let 
f G B(F), and a < b. Lf f is F -continuous on F PI [a, b], then f is F a -integrable on 
[a,b]. 

Proof: If Sp(a) = Sp(b), then the F Q -integral is zero, and the result is obvious. 

Now consider the case where Sp(a) ^ S F (b). The function / is uniformly F- 
continuous on F n [a, b] as F n [a, b] is compact. Thus, given e > 0, there is a 5 > 
such that 

x, y G F n [a, b] and \y - x\ < S =>■ \f(y) - f(x)\ < ^ ^ g ^ a y (34) 

Let P be a subdivision such that \P\ < 6. Then it can be seen that U a [f,F,P] < 
L a [/, F, P] + e which completes the proof in view of lemma OH • 

The following property of ^"-integral is expected from any fair definition of an 
integral: 

Theorem 40 Let a < b and f be an F a -integrable function on [a,b]. Let c € (a,b). 
Then, f is F a -integrable on [a,c] and [c, b]. Further, 

f(x) d a F x = f f(x) d a F x + t f{x) d a F x (35) 



This can be proved in a manner analogous to Riemann integral. 
The linearity of F"-integral follows from the definition: 

Theorem 41 (i) Lf f is F a -integrable on [a,b], and A is any real number, then 



Xfix)d F x = X f(x)d F x. 

J a 

(ii) If f and g are F a -integrable functions on [a,b], then 

ifix)+gix))d F x= / fix)d a F x+ I gix)d F x. 



The following lemma states an obvious property: 

Lemma 42 If f and g are F a -integrable over [a,b], and fix) > gix) for all 
x G F (~1 [a, b], then 

fix)d a F x> f gix)d a F x 

J a 

Definition 43 If f is F a -integrable on [a,b], a < b, then 

,6 

fix)d F x = - / fix)d F x. 



A particularly simple but important example, as realized in |16| . is the ^"-integral 
of the characteristic function xf of the set F: 
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Lemma 44 If Xf{ x ) is the characteristic function of F cR, then 

fb 

X f{x) d F x = S F (b) - Sf(a) 



Proof: For a closed interval / C [a, b], 

M[xF,F,I]=m\xF,F,I} = l ifFn/^0 

= otherwise 

so that M[xf, F, I] is zero for a closed interval / = [c, d] only when Sp(d) — Sp(c) = 0. 
Thus 

U a [xf ,F,P] = L a {xf ,F,P] = Sp(b) — Sp(a) 
for any subdivision P of [a, b]. • 

As a further example, in | Appendix A| we calculate the C Q -integral of the function 
f{x) — xxc{x) where C is the middle | Cantor set, and a — ln(2)/ln(3) is its 
7-dimcnsion. The integral is given by equations (|A.12|) and (|A.13I) . 

8. F Q -DifFerentiation 

Like the first order derivative, the ^"-derivative is a limit of a quotient. But here the 
limit is f-limit, and the denominator is the difference in the values of the staircase 
function Sp at two points. Moreover, intuitively speaking, F is typically the set of 
change of the function, and a is typically the 7-dimension of F. 

Definition 45 If F is an a-perfect set then the F a -derivative of f at x is 

Uhm c m r{ (X ), itxtF 
V F (f(x)) - <^ »-* S F (y) - S F (x) (36) 

I otherwise 

if the limit exists. 

Note that lemma I2TH tells us that if x € F, then we would find such points y which are 
arbitrarily close to x at least on one side of x so that the denominator in the definition 
is not zero and the RHS in (|36|l makes sense. 

We now state a necessary condition for the above limit to exist. 

Theorem 46 IfT> F {f{x)) exists for all x in (a, b), then f(x) is F -continuous in (a, b). 

The proof is straightforward. 

The linearity of the F Q -derivative is an immediate consequence of the 
definition 23 Thus: 

Theorem 47 (i) Let f be a function on [a,b]. If V F (f(x)) exists for all x G [a, b], 
then T> F (\f(x)) exists and 

V F (Xf(xj) = XDp{f{x)). 

(ii) Let f and g be functions on [a,b]. IfD F {f{x)) andT> F {g(x)) exist for all x £ [a, b], 
then T> F (f(x) + g(x)) exists and 

VfV(?) + g(x)) - V F (f(x)) + V F {g(x)). 



Calculus on fractal subsets of real line - I: formulation 21 

Now we calculate the derivative for two elementary functions. The first does not 
need a proof: 

Lemma 48 The F a -derivative of a constant function f : R — > R, f{x) = k £ R is 
zero, i. e. 

V F {f) = 0. 

This result is to be contrasted with the classical fractional derivative (Riemann- 
Liouville, and others) of a constant which is not zero in general [291 131)1 13"T1 132) . 

Lemma 49 The derivative of the integral staircase itself is the characteristic function 
XF ofF: 

V f (S%(x))= X f(x). 

Proof: If x i F, V F {S^{x)) = 0. 
If x £ F, then 

Vf(S${x)) = F-\to fjf\~ fjf\ = 1 



This lemma together with lemma 2H can be viewed as the special cases of the 
fundamental theorems of calculus (section 0) involving S F and its derivative xf- 
An analogue of Rolle's theorem is: 

Theorem 50 Let f : R — > R be a continuous function such that Sch/ C F where F 
is a-perfect, T> F (f(x)) is defined for all x £ [a, b], and /(ft) = f(b) = 0. Then there 
is a point c £ F n [a, b] such that D?(/(c)) > and a point d £ F n [a, 6] such that 
Vf(f(d))<0. 

Proof: If / is zero throughout [a, b], then V F (f(x)) = for all x £ [a, b] and the result 
follows in this case. 

If f(y) > for some y £ (a,b), then as / is continuous, there exists an open 
interval (c, d) C (a,b) containing y such that f{z) > for any z £ (c, d). Let (co,do) 
be largest such interval. Then /(cq) = /(do) = 0. The point cq £ Sch/ C F as / is 
positive on the right of cq. So, 

V F (c ) = F-Mm 9 ff Z j-ffi\ >0. 
z^co b%(z) - <S£(c ) 

Similarly do £ F and 

with the same considerations. The points Co and do can be identified as points c and 
d in the statement of the theorem. 

If there are no points y such that f(y) > and neither is the function zero 
throughout, then we can choose a point y such that f(y) < and proceed in a similar 



Remark: The following example shows that the analogue of Rolle's theorem can not 
be made more strict which would have implied existence of a point c £ F such that 
V F (f(c)) = 0. Let C be the middle \ Cantor set. Define 

f(x) = Sc(x) < x < 0.5 

= l-5g(x) 0.5<x<l. 
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This function satisfies /(0) = /(f) = 0. Further, it is continuous in the interval [0, f]. 
Its set of change is C. The C Q -derivative is given by 

VcU(*))=Xc{x) 0<x<0.5 
= -xc(x) 0.5 < x < 1. 

Thus, x G C => Pg(/(a;)) = ±1^0. 

In general, it can be said that the "fragmented nature" of the fractal F does not 
allow us to make the analogue of Rolle's theorem as strict as its original version. 

Now we state the analogue of the law of the mean. 

Corollary 51 Let f : R — > R be a continuous function such that its set of change 
is contained in an a-perfect set F C R, V F {f{x)) exists at all points x G [a, b] and 
Sp (b) ^ Sp (a) . Then there exists a point c G F such that 

f(b) - /(«) 



mm) > 

and a point d G F such that 
V a F (f(d)) < 

Proof: Let 



S%(b)-S F (a) 

f(b) - f(a) 
S F (b) S F (a) 



g(x) = (f(x) - /(a)) - ^M_|^L(5^(x) - S£(a)) 

so that the difference between / and g is a constant plus a multiple of Sp (x) — Sp (a) . 
Now g(a) = g(b) = so that we can use the last theorem to say that there exists a 
point c G F such that V F (g(c)) > and a point de F such that V F (g(d)) < 0. This 
implies that 

for some c G F, and 

VAf{d)) Sf(b)-Sf(a)-° 
for some deF, which lead to the required relations. • 

We had seen earlier that the ^"-derivative of a constant f(x) = k is zero. Now 
we see that these are the only functions whose ^-derivatives are zero: 

Corollary 52 Let f : R — > R be a continuous function such that Sch(/) C F and 
V F (f(x)) — for all x G [a, b]. Then f(x) = k where k is a constant on [a, b}. 

Proof: Suppose, if possible, that the function is not a constant. Then there exist y 
and z, y < z, such that f(y) ^ f(z). This implies either f(y) < f(z) or f(y) > f(z). 
Case 1. f{y) < f(z). Then there exists a c G F PI (y, z) such that 

Case 2. f(y) > f(z). Then there exists a d G F n (y, z) such that 
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In both the cases we have found a point where the derivative is not zero which 
contradicts our assumption. • 

Remark: Again due to the "fragmented nature" of the fractal F, the F a - 
differentiability of / is not sufficient to guarantee the result. Further, The additional 
conditions that / be a continuous function and Sch(/) C F are necessary also in the 
second fundamental theorem [SSJ which relies on the last corollary, and the integration 
by parts rule (theorem [SHI) which depends on theorem 1551 
The F"-derivative satisfies the analogue of Leibniz rule: 

Theorem 53 If the functions u : R — > R and v : R — > R are F a - differ entiable, then 
fix) — u(x)v[x) is F a - differentiate, and 

V F (u(x)v(x)) = V F (u(x))v(x) + u(x)V%(v{x)). (37) 

The proof is straightforward. 



9. Fundamental theorems of F"-calculus 



This section relates the F a -integration and ^"-differentiation as "inverse processes" 
of each other. 

The first fundamental theorem says that the ^"-derivative is the inverse of 
indefinite F"-integral. 

Theorem 54 Let F C R be an a-perfect set. If f € B(F) is an F '-continuous 
function on F H [a, b] and 



g(x)= / f(y)d F y 

J a 

for all x G [a, b], then 

V%(g(x)) - f(x) XF (x). 

Proof: If x $ F, then T>p(g(x)) — by definition. 

For x S F, if there are points in F arbitrarily close to x on both sides of x, then 
we have to consider both the following cases: 

(i) The set F n (x, z) is never empty for z > x and 



g(z)-g(x)= / f(y)d a F y. 

J X 

(ii) The set F D (z, x) is never empty for z < x and 
g(x)-g(z)= [ X f(y)d a F y. 



Otherwise we have to consider only one of the cases which is applicable. We consider 
the first one; the second can be treated similarly. 

In the first case, F n (x, z) is not empty for any z > x. Taking the F-limit as 

z — » x, we get 

f z f(v) d a v 

Now, 

m[f,F,[x,z]] f XF{y)d F y<( f(y)d F y<M[f,F,[x,z}][ XF {y)d F y, 

J X J X J X 
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and 

f XF{y)d F y = [Sp{z) - S a F {x)) 

J X 

so that 

f z f(y) d a y 

m[L F ' [X ' Z]] - S«(z ) - S°(x) ~ M[L Fl [x ' Z]l m 
As / is F-continuous, 

F- lim m[f, F, [x, z]] = F- lim M[f, F, [x, z}] = f(x) (40) 

z — >x z — >x 

From |J5SJ|, (|39|l and (|40|) . we get the required result. • 

The second fundamental theorem says that the ^"-integral as a function of upper 
limit is the inverse of ^"-derivative except for an additive constant. 

Theorem 55 Let f : R — > R be a continuous, F a -differentiable function such that 
Sch(/) is contained in an a-perfect set F and h : R — + R be F-continuous, such that 

h(x) X F(x)=V F (f(x)). 



The 



Proof: If 



then V F (g(x)) = h{x)XF{x) by the last theorem. Therefore V F (g(x) — f(x)) = 
for all x 6 [a, b]. Now corollary 1521 implies that g(x) — f(x) — k, a constant, or 
g{x) = f{x) + k. Thus, 

fib) - f(a) = g(b) - g(a) = g(b) = [ h(x) d F x 



b 

h(x)d F x = f(b)-f(a). 



g(x) — / h(x) d F x 



which proves the theorem. 



The following theorem states that the ^"-integration can be performed by parts, 
and can be proved by using fundamental theorem 1551 and Leibniz rule ftheorem 153(1 : 

Theorem 56 Let the functions w:R— >R, u:R— >R be such that 

(i) u is continuous on [a, b] and Sch(u) C F, 

(ii) T> F {u) exists and is F-continuous on [a,b], 
(Hi) v is F-continuous on [a, b]. 

Then, 



h 



lb rb 



u(x) / v(x')d F x 



V F {u{x)) / v(x') d F x' d F x. (41) 



u(x)v(x) d F x 

The proof is straightforward and omitted. 

In | Appendix B| we discuss examples of repeated F"-derivatives and ^-integrals. 
There we also calculate F"-derivatives and ^"-integrals of powers (S F {x)) n . 

Now that the analogies between F Q -calculus and ordinary calculus have become 
clear, we summarise some of them in section [Appendix C| for a quick reference. 
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10. Examples and applications of F Q -Differential equations 

In this section we briefly touch a couple of examples of F Q -differential equations. The 
^^-differential equations is the main topic of a subsequent work [3"§| . 

Firstly we revisit the local fractional diffusion equation proposed in and also 
discussed partly in |17j . This equation is of the form 

V% tt (W(x,t)) = *?&^W(x,t). (42) 

where the density W is defined as a function of two arguments (x, t) 6 R x R and 
with a slight change of notation Dp t denotes the partial ^"-derivative with respect 
to time t, xf being the characteristic function of F. (This equation may be compared 
with ordinary diffusion equation ^gf^ = D-J^W(x, t).) The Riemann integral like 
prescription given in |16j had enabled one to construct a new exact solution. This 
solution is 

W(x,t) = rexpf— ^— V W(x,0) = S(x). (43) 

V ' {2irS%(t))i P \2Sp(t)J K ' ' y ' y ' 

This can be recognized as a subdiffusive solution, since Sp is known to be bounded 
by kt a , k constant, in simple cases including Cantor sets. 

An important observation at this stage is that: equations like (14211 are examples 
of fractal-time evolution processes. 



Motion in a fractally distributed medium 

As a second example, we consider one dimensional motion of a particle undergoing 
friction. First we recall the equation of motion in a continuous (i. e. nonfractal) 
medium. If the frictional force is proportional to the velocity, the equation of motion 
can be written as 

^ = -Hx)v (44) 

where k(x), the coefficient of friction, may be dependent on the particle position x. 

Equation (|44|l can be reexpressed by considering velocity v as a function of 
position x. The equation can be written as 

dv dx 

— = —kixjv. 

dx dt y ' 

Identifying dx/dt = v and assuming d^O, the equation becomes 

^ = ~Hx) (45) 
da; 

which is readily solved by integrating k(x) if k(x), which models the frictional medium, 
is smooth. 

If the underlying medium is a fractal, then (|45(l is inadequate to model the motion. 
Instead we propose the ^"-differential equation of the form 

Vp{v(x)) = -k(x) (46) 

for this scenario. Here, the set F is the support of k(x) which describes the underlying 
fractal medium, and a is the 7-dimension of F . (If F is not a-perfect, then the set 
Sch(S^) can be chosen instead.) The function k(x) may be called fractional coefficient 
of friction due to its physical dimensions. 
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The solution of il46|) is easily seen to be 
v(x) — vq — / k(x r )dpx / 



(47) 



where ?;o and xq are the initial velocity and position respectively. In a simple case 
where k{x) is uniform on the fractal i. e. k(x) = kxf(%) where k is a constant, (|47H 
reduces to 



In the extreme cases we obtain back the classical behaviour: (i) If F is empty 
(frictionless case), then v(x) — vq; (ii) If F = R (uniform medium) then v(x) = 

Vq - k(x - X ). 

The time dependence of x is given by 



where t(x) is the time required to reach the position x. 
11. Concluding remarks 

In this paper we have developed a calculus on fractal subsets of the real line. 
This developement involved the identification of the special role played by staircase 
functions associated with fractal sets, which may be compared with the role of 
independent variable itself in ordinary calculus. In particular, ^"-integrals and F a - 
derivatives (of order a, < a < 1) are defined using staircase functions for sets F 
of dimension a. In contrast with the classical fractional calculus, the notions of F a - 
derivatives and ^"-integrals are specifically tailored for fractals of dimension a and 
thus provide suitable operators on fractals. Further, they reduce to ordinary derivative 
and Riemann integral respectively, when F — R and a = 1. 

Much of the developement of the ^"-calculus is carried in analogy with the 
ordinary calculus. Several results and techniques of ordinary calculus, including the 
Leibniz rule, the fundamental theorems of calculus, the technique of integration by 
parts etc. have analogues in this calculus. Specifically we have adopted Riemann 
approach for ^"-integrations. This approach can possibly be generalised using 
Kurzweil-Hcnstock integration schemes |40l 141) . Work is in progress in this direction. 

In the process of the developement of the ^"-integrals we have introduced a- 
mass or mass function associated with a fractal subset F of the real line. This lead 
us to introduce the 7-dimension in section [3] This dimension is finer than the box- 
dimension. Though it is not as fine as the Hausdorff dimension, it is specific to 
the developement of calculus here and we expect it to be associated naturally with 
algorithms and numerical schemes based on the present calculus. 

We have also discussed simple models based on ^"-differential equations. 
The solutions of ^"-differential equations naturally involve staircase-like functions. 
Staircase functions such as the Lebesgue Cantor staircase function are known to 
be bounded by sublinear power laws. Also, they "change" or "evolve" only on a 
fractal set. Thus, this framework may be useful in modelling many cases of sublinear 
behaviour, fractal time evolution, fields due to fractal charge distributions, etc. The 
^"-differentiability may be used to classify singular probability distribution functions. 

Continuous-time dynamical systems are associated with ordinary differential 
equations, and discrete-time dynamical systems are associated with maps/ 



v(x) =v - K(Sp(x) - Sp(x ))- 
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diffeomorphisms. But as realised in the dynamical systems associated with F a - 
differential equations would be those evolving on fractal subsets of time-axis. It would 
also be of great interest to investigate correspondences between ordinary differential 
equations and F a -differential equations. These are explored in a companion paper |39 |. 

There are many obvious directions in which considerations of this paper should be 
extended. Some of them are mentioned above. Other important directions would be 
extensions to multivariable case, developement of differential equations and variational 
principles to mention a few. Work is in progress in these directions. 
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Appendix A. C"-Integrating f(x) = xxc(x) 
As an example of ^"-integration, here we calculate 



g(y) = / xxc(x)dcx 



xd^x 



(A.l) 



where C is the middle | Cantor set, and a — ln(2)/ln(3) is its 7-dimension. The 
function f(x) — xxc{x) is C-continuous on [0, 1], hence it is C Q -integrable. 

The set P n = {xi = i/n : < i < n} is a subdivision of [0,1]. For any 
component [xt,x l+ i] of P n , Xi < m[f, F, [xi, and Xi+i > M[f, F, [n, Xi+i]] if 

F n [xi,Xi+i] ^ 0. Therefore, 



7(1) = km £ 



i=0 



QCX 



5g 



<L a [f.F,P n 



and 



^(1) = lim J2 



s 1 ; 



i + l 



S° r 



>U a [f, F,P n }. 



Further, it can be seen that 



lim [.g(l) - g(l)} = 0. 



Thus, 3(1) can be calculated using the limit 

'i + l 



9 



(1)= lim £ -L 



5? 



(A.2) 



Similarly for integers m > 0, 



9 



— ) = lim y 



3 m n 



1 



3 m n 



3 m n 



(A.3) 



Using the self-similarity of C and scaling of Sq, we see from (|A.2|) and (|A.3|I that 
1 



9 



1 1 

3m(l+a) 9{ 1) — ^77 .9(1)- 



(A.4) 
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We make use of the ternary representation of numbers which simplifies many 
calculations involving the Cantor set. Any number y S [0, 1] can be represented by 
the series 

8=1 

where U(y) — 0, 1 or 2 is the ith ternary digit of y after ternary point. The number y 
belongs to C if and only if y has a representation of the form (|A.5J) where (y) = or 2 
for all i. 

An approximation of y G [0, 1] by a finite number of digits is denoted by 

T o (y) = and T n {y) = £ (A.6) 

i=l 

The sequence {T Tl (y)}^°_Q is a monotonically (but not strictly) increasing sequence 
whose limit is y. Hence we can write 

oo pTi(y) oo 

g(y) = J2 *x(*)d£x = 5>(y) (A.7) 

where 

h(y)= / x X {x)A a c x. (A.8) 

The quantities Ii(y) can be calculated using the self-similarity of C, the scaling and 
translation properties of S*g (theorem I14f) . and (|A.4|I . Let y G [0, 1] and let n be any 
integer such that i < n ti(y) = or 2. Then i < n =4* Ti(y) G C. For calculating 
/„ (y) , there are three cases corresponding to three possible values of t n (y) : 
Case t n {y) = 0: Here, T n - X {y) = T n {y) and I n {y) = 0. 
Case t n (y) = 1: In this case, 

1 00 2 
T n (|/)-T n _i (|/) = — = # 

i—n+l 

so that there is another sequence {ti(T n (y))} which does not contain the digit 1, 
hence T n (y) G C. The set [T n - 1 (y),T n (y)} n C can be written as 

{z : i <n i;(z) = ^(y); t n (z) = 0; i > n =>• ^(z) = or 2} 

Therefore it is a scaled down version of C by a factor 1/3™ and translated by 
T n -i(y). Hence writing x = T„_i(y) + (x - T n -i(y)), we get 

I n (y)=T n - 1 {y) / Xc{x)d a c x+ {x-T^iy^xcix)^ 

rl/S" ,1/3™ 

= T n -i{y) / xcW4 ;I + / a^Xc(a;)dg.a; 
Jo Jo 

= + 3 n(i+a) g( x ) 

_ T»_ x (y) g(l) 

where T denotes r(a + 1) for convenience. If y = T n (y) then y 6 C. But if 
2/ > T n (y), then as i n (y) = 1 and £; 7^ for some i > n, therefore y ^ C. Thus 
the half open interval (T„(y), y] does not intersect C implying that Ife(y) = for 
all k > n. 
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Case t n (y) = 2: Here, T n (y) clearly belongs to C. If D is the set 

D = {z : i < n => U(z) = U(y); t n (z) — 0; i > n => ti(z) = or 2} 

then D is a scaled down version of C by a factor 1/3™, D c [T„_i(y), T n (y)], and 
more specifically, [T n -i(y),T n (y)] n C = D U {T„(j/)}. Therefore by arguments 
similar to the case t n (y) = 1, 

/„(,) = + (A.9) 

But unlike the case £ n (y) = 1, there is a possibility that CCi(T n (y),y] is nonempty 
so that Ik{y) need not be zero for all k > n. 

Summarizing, 

ift„(j/) = 0or 

/.,(/,)= < t,(y) = 1 for some i < n (A1Q) 

i^m + m. otherwise . 

p 2 « Qn 

This description requires the value of <?(1). It can be found out by putting y = 1 in 
{OOjl. If y = 1, then t t {y) = 2 for all i. Also, 

n „ 

i=i d 

Therefore, 

z n( i) = 1 - 3 " ( "" 1) + ^. 

" v y T2" 6™ 

Substituting this in (|A.10|) and solving (| A. T|> for <?(1), we get 

5(1) = ^- (A.ll) 



Thus, 



where 



g(y)= x Xc (x)d a F x = J2in(y) (A.12) 

^0 n=l 



1 



r(a + l) 



r„-i(y) 



2 • 6™ 



if i„(y) = or 

ii(y) = 1 for some i < n 

otherwise 



(A.13) 



and T n (y) are given by equations (| A. 5|> and i|A.6|) . 



Appendix B. Regarding repeated ^"-integration and F Q -derivative 



F a -derivative 

Many dynamical systems are modelled by differential equations involving second and 
higher order derivatives i. e. derivative operator applied repeatedly. The successive 
operation of the T> F operator is also possible and gives meaningful results. As an 



Calculus on fractal subsets of real line - I: formulation 30 

example, let us ^"-differentiate the function g(x) = (Sp(x)) 2 twice, where F C R is 
an a-perfect set. By definition of the derivative, 

x i F => V F g(x) = 0. (B.l) 

If x G F, then 

P^-F-hm 5^ _ 

= 2S£(x). (B.2) 
Equations i|B.l|) and (|B.2|1 can be combined to give 

V a F (S£(x)) 2 =2S%(x) X f(x) (B.3) 

where \f is the characteristic function of i* 1 . As a side remark, it is easy to generalize 
this to 

V F ((S^x)) n ) = niSftx^XFix) (B.4) 

for any integer n > 0. 

Now we take the second F a -derivative of g. As far as the operator T> F is 
concerned, the values of the function outside F make no difference because of the 
i^-limit in its definition. Thus, 

(V f ) 2 (S^(x)) 2 ^V f (2S^x) X f(x)) 
= 2V F S%{x) 

= 2 X f(x) (B.5) 

Where the last step follows from lemma RUl and linearity (theorem 1470 . 

Apart from the 7-dimension of F, the order a also has another significance. This 
will be clear from the following example. If C is the Cantor set, then it is known 42 
that Sq{x) is bounded by the power a of x from below and above: 

ax a < 5g(x) < bx a (B.6) 

where a = ln(2)/ ln(3) is the 7-dimension of C . Hence the function g defined above is 
bounded by power 2a of x: 

ax 2a < g{x) = (S%(x)) 2 < bx 2a (B.7) 



so that 



and 



x £ F 2ax a < V%{g{x)) < 2bx a (B. 



xeF^(V%) 2 (g(x))=2. (B.9) 
This example demonstrates that F Q -differentiation reduces the power of bounds by a. 
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Table CI. A few analogies between F^-calculus and ordinary calculus. 



Ordinary calculus 



.F^-calculus 



R 

limit 

Continuity 
x n Ax = 



n + 1 

— x n =«/-' 
dx 

Leibniz rule 
Fundamental theorems 
Integration by parts 



An a-perfect set F 

F-limit 

_F-continuity 

' '{Sf(x)) n A a F x = 

'0 



1 



n+1 



(S F (y)T 



-D%((S%(x)) n )=n(S%(x)r- 1 X F(x) 
Theorem 1531 
Theorems 1541 and 1551 
Theorem 1561 



F a -integration 

The ^"-integration can also be carried out in succession. Let F be an a-perfect set. 
It is already shown that F"-integration of Xf(%) is S F (x): 



X f(x) d a F x = S F (x') 



(B.10) 



where for simplicity we have taken S F (a) = which is consistent with the definitional 
of the staircase function. Now, 



91 (*') 



S F (x)d F x = ~(S F (x')) 2 



where we have used (|B.3|I and the fundamental theorem [331 Again, it is easy to 
generalise this to 



(S F {x)) n d F x 



(s F (x r )r 1 



n + 1 



(B.ll) 



using (|B.4I 



Appendix C. A few analogies between F"-calculus and ordinary calculus 

The F Q -calculus can be thought of as a generalization of ordinary calculus with 
Riemann approach. Table shows a few analogies between various quantities. 
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